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ABSTRACT 

Context. A wealth of observations of CO in absorption in diffuse clouds has accumulated in the past decade at uv and mm-wavelengths 
Aims. Our aims are threefold: a) To compare the uv and mm-wave results; b) to interpret "CO and '-CO abundances in terms of the physical 
processes which separately and jointly determine them; c) to interpret observed J=l-0 rotational excitation and line brightness in terms of 
ambient gas properties. 

Methods. A simple phenomenological model of CO formation as the immediate descendant of quiescently-recombining HCO+ is used to study 
the accumulation, fractionation and rotational excitation of CO in more explicit and detailed models of H2-bearing diffuse/H I clouds 
Results. The variation of N(CO) with N(H2) is explained by quiescent recombination of a steady fraction n(HCO^)/n(H2) = 2 x 10"'. Observed 
N('^CO))/N('^^CO) ratios generally do not require a special chemistry but result from competing processes and do not provide much insight into 
the local gas properties, especially the temperature. J=l-0 CO line brightnesses directly represent N(CO), not N(H2), so the CO-H2 conversion 
factor varies widely; it attains typical values at N('^CO) S 10'^ cm"^. Models of CO rotational excitation account for the line brightnesses 
and CO-H2Conversion factors but readily reproduce the observed excitation temperatures and optical depths of the rotational transitions only if 
excitation by H-atoms is weak - as seems to be the case for the very most recent calculations of these excitation rates. 

Conclusions. Mm-wave and uv results generally agree well but the former show somewhat more enhancement of '^C in '^CO. In any case, 
fractionation may seriously bias '^C/'^C ratios measured in CO and other co-spatial molecules. Complete C— »CO conversion must occur over a 
very narrow range of Ay and N(H2) just beyond the diffuse regime. For N(H2) < 7x10" cm"^ the character of the chemistry changes inasmuch 
as CH is generally undetected while CO suffers no such break. 
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1. Introduction. 

Except for hydrogen, carbon monoxide is the most impor- 
tant and widely observed molecule in the interstellar medium 
(ISM). The 7-8 decade span in column density over which CO 
is directly observed, from N(CO) = lO'^ cm"^ in uv absorp- 



gas-phase ion-m olecule chemistry of lo w -dens it y media like 



diffuse clouds jVan Dishoeck & Blackl 11988c iWarin et al 
1996). The relative abundance of CO with respect to H2 
varies widely in this regime (see Fig. 1; with much scatter, 
approximately as N(CO) oc N(H2)^ over the range X(CO) - 
N(CO)/N(H2) ~ 3 X 10-** - 3 X 10-^) and the relative abun. 



tion in the diffuse interstellar medium to N(CO) > 10^^ cm ^ ■ ei2r^r\ a i3r^r\ * 1 ^ * * 

^ , dances of CO and CO are strongly affected by fractionation 



in mm and sub-mm emission from dense and giant molecular 
clouds, is exceeded only by that of H2 itself. The ubiquity of 
CO has encouraged the use of mm-wave CO emission as a pos- 
sible tracer of molecular hydrogen even into such extreme en- 



jWatson et all Il976t ISmith & Adamsi Il980l) such that 20 ^ 



vironments as high velocity clouds (iDessauges-Zavadsky et al 
2007h . 

Interpreting observations of CO in diffuse gas 



N('2cO)/N('^CO) ^ 170 (see the references cited in Sect. 2) 



Despite the comparatively small CO abundances in diffuse 
clouds, A 2.6mm J=l-0 rotational emission is often apprecia- 
ble. Typically it is seen that Tb ~ 1 -5 K, Wco~ 1 -5 K km s"^ 

-2 



dSnow & McCallL 12006.) over the lower half of its range, peak or integrated br ightness for N(CO) = lO'^ - lO"" cm 



at N(CO) ^ 10 cm , has been particularly challenging. 
The fraction of free gas-phase carbon in CO is small, a 
few percent or less when, locally, Av< 1 mag, but it is still 
30-50 times larger than can be explained by the quiescent 
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( Liszt & Lucas , 1998|) but peak brightnes ses as high as 10-13 K 
have been observed Liszt & Lucasl(ll994l) . Moreover, the J=l-0 
rotational transition may have appreciable optical depth in dif- 
fuse gas because its excitation is quite weak. Surveys of '^CO 
alone may be hard-pressed to distinguish between dark and dif- 
fuse gas, especially at higher galactic latititude or larger galac- 
tocentric radii owing to the broader distribution of diffuse gas. 
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To elucidate the properties of CO in the diffuse regime, 
we discuss here a weahh of observational material at uv and 
mm-wavlengths which has accumulated (much of it very re- 
cently) over the past decade. The plan of this work is as 
follows. Sect ion 2 gathers the previously-published observa 



tional results (Liszt & Lucas'.'1998':'S onnentrucker et al.L 12007 



[Burgh et al., 2007; Sheffer et al., 2003) which form the basis 
of the present discussion. Section 3 displays and discusses the 
run of observed values of the CO and H2 column densities to 
demonstrate that there is at least a phenomenological basis for 
understanding the abundance of CO in diffuse gas, in order to 
show that there is some knowledge of the microscopic CO for- 
mation rate. This rate and those of the various other physical 
processes which account for the abundances of '^CO and '^CO 
are set out in detail in Sect. 4. Section 5 discusses the observed 
abundances of '^CO and '^CO and their fractionation, and Sect. 
6 discusses the rotational excitation and brightness of the J=l- 
rotational transition. Section 7 is a brief discussion and sum- 
mary of outstanding concerns. 



2. Observational material 

With one minor exception (see Sect. 6) , the discussion here 
rehes on previously-published results in the uv and mm- 
wavelength regimes, as we now discuss. 



2.1. uv and optical absorption 

Many determ inations of N(CO) and N(H 2) have recent l y been 
pub lished by ISonnentnicker et al.l (120071) . iBurgh et alj (120071) 
and lSheffer et al.l (l2007^. who provide measurements of N(H2), 
N('^CO) and N(^^CO), along with such additional physically 
interesting quantities as the J=l-0 rotational temperatures of H 2 
and/or carbon monoxide and, from So nnentrucker et al. ( 2007 ), 



column densities and related quantities for such species as CH 
and C2. For the lines of sight where there is overlap, agree- 
ment is generally excellent for the column densities of CO 
and H2 and only slightly worse for the CO excitation temper- 
ature, as shown below in Figs. 5 and 6. In cases of overlap, 
we chose values from the reference with the smaller quoted 
errors if the datasets were equally comprehensive in that dkec- 
tion. However, we also chose not to mix values for N('^CO) 
and N('^CO) from different references along any given line of 
sight, because the systematic errors could be different. 

2.2. MM-wave absorption 

Observations of carbon monoxide in a bsorption towa r d mm - 
wave continuum sources were given by L iszt & LucasI (119981) . 
along with rotational excitation temperatures and isotope ra- 
tios, etc.). N(H2) is not known directly in these measure- 
ments, which a ssume instea d that N(H2) - N(HC0+)/2 x 10"^ 



dLiszt & Lucaslll 996: Lucas & Lisztlll996i,,2000.) . The cuiTent 



discussion may be regarded as a consistency check on this as- 
sumption. 

In comparing the radio and optical lines of sight, it should 
be remembered that the former use extragalactic background 




^Gredel 
owelty 



N(H2) [cm-2] 



Fig. 1. Top: variat ion of CO and H2 co l umn d ensit ies observed 
in uv ( rectangles; Sonnentrucker et aTl (2007h and Burgh et al 



(I2OO7 I)) and mm- wave (triangles: iLiszt & LucasI (Il998l) ) CO 



absorption. Lines of sight represented in Fig. 5 are outlined and 
shown in green. For the mm-wave data N(H2) = N(HCO^)/2 x 
10"^. The curves represent models of CO formation via recom- 
bination of HCO^ and are labelled by their density n(H) ^ n(H 
I) + 2n(H2) (see Sect. 4). The lowest curve shows the result 
of halving the HCO^ abundance at n(H) =16 cm"^. Bottom: 
Variati on of CH and H? column de nsities from the summary ta- 
bles of 'Sonnentrucker et al] (l20C)7l) ar id a few high-N(CH ) dat- 
apoints from Gredel et a]T(ll993ir and lWelty et alJ (l2003h . 



sources, penetrate the entire galactic layer, and refer to individ- 
ual, well-resolved kinematic components. The optical/wv lines 
of sight stop within the Galaxy and are sums over undifferen- 
tiated - though not necessarily blended - features for species 
studied in the uv, that is for CO and H2. 
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3. Variation of molecular abundance with H2 and 
the source function for carbon monoxide 

The abundances of '^CO and '^CO are the result of several pro- 
cesses; chemical formation and destruction, carbon isotope ex- 
change and (selective) photodissociation. Interpreting the ob- 
servations requires an understanding of the competing influ- 
ences of all such processes, but the overall guage is really set by 
a comparison of the photodissociation rate (which is known in 
free space and calculable within a model) and the direct chem- 
ical formation rate. Once the latter is specified, other physical 
properties follow directly and it is even possible to model the 
internal rotational excitation and predicted mm-wave bright- 
ness. Here we outline an empirical approach to estimating the 
rate of CO formation. 

3.1. C0andH2 

In Fig. 1 at top we show the observed run of N('^CO) with 
N(H2). The rise of N(CO) with increasing N(H2) is very steep, 
approximately N(CO) oc N(H2)^, with scatter amounting to two 
orders of magnitude at fixed N(H2). At each N(H2) the radio- 
derived N(CO) is comparatively but not inordinately large; this 
could be an indication that X(HCO^), taken as 2 x 10 has 
been slightly overestimated. The large subset of wv-absorption 
datapoints which have CO excitation temperature measure- 
ments and which therefore appear also in Fig. 5 have been 
flagged and colored green in Fig. 1, in order to highlight the 
fact that most of them require model densities of at least n(H) 
= 64 cm"-' (see Sect. 6). 

What to make of the variation of N(CO) with N(H2)? 
Models of CO formation employing thermal processes in qui- 
escent diff'use gas have largely been unable to reproduce ob- 
served values of X(CO) except perhaps at rather higher densi- 
ties than are otherw i se infe rred for diffuse gas; see Fig. 16 of 
Sonnentrucker et al.l (l2007h . However, purely phenomenologi- 
cally, it is known that thermal gas-phase electron recombina- 
tion of the observed amount of HCO+, X(HCO+) ~ 2 x lO '', 
does suffice to reproduce the observed X(CO) at modest den- 
sities. This is shown by the curves at the top in F ig. 1, where, 
updating the calculations of iLiszt & Lucas poooV we plot the 
predicted N(CO) for diffuse cloud models in which an artifi- 
cially steady relative abundance X(HCO^) = 2 x 10 is al- 
lowed to recombine with free electrons at the thermal rate; note 
that the mm-wave observations of N(CO) and N(HCO^) have 
been placed in the plane of the Figure by relying on the same 
X(HCO^) which in the models explains the CO abundance. 

The models are for small uniform density spheres im- 
mersed in the mean galactic radiation fields, within which the 
thermal balance, ionization equilibrium and H2 and_CO accu- 
mulation problems are solved sel f-consistently ( Liszti l2007b 
using the self-shielding factors of iLee et al.l (119961) . For each 
total density n(H) = 16, 32, 64 ... 256 cm"-', a series of models 
of increasing N(H) was calculated, and the CO and H2 column 
densities across the central line of sight were calculated to form 
the curves for each n(H). 

The newer calculations differ from the older ones in a va- 
riety of small ways, for instance the rates for O I and C I ex- 



citation by hydroge n atoms were very recently recalculated by 
Abrahamsson et al. (EootI) . VIost important is the recognition, 
following observation of H^^, and accounting for the neutral- 
ization rate on small grains/PAH, that that the low-level hard 
ionization rate of hydrogen, presu mably due to cosmic rays , 
is apparently larger in diffiise gas (iMcCall et al.L l2002t iLiszi 
2003). This has the effect of increasing the densities of H^ and 
He^ and the overall electron fraction (which, however, is never 
more than about twice the free carbon abundance). Higher elec- 
tron fractions lead to more rapid formation of CO (if CO forms 
by HCO^ recombination) while the increased density of He^ 
plays an elevated role in the destruction of carbon monoxide. 
Indeed, destruction by He^ is the dominant chemical mecha- 
nism for '^CO destruction. The various processes responsible 
for destroying CO are discussed in Sect. 4. 

In the Figure, scatter in N(CO) at fixed N(H2) is implicitly 
attributed to variations in density, which in the models ranges 



between n(H) =16 cm and 256 cm . However, the very 
lowest curve shows the result of halving X(HCO^) in the model 
and additional observational scatter , perhaps substantial, could 
also result from variations in geometry and ambient illumi- 
nation. Consideration of fractionation strongly suggests vari- 
ations of a factor two or more in the photodissociation rate (see 
Sect. 5 and Fig. 4). 

In the remainder of this work we will assume that the direct 
CO formation rate is adequately specified by the dissociative 
recombination of HCO+ at X(HCO+) = 2x10"'' and we will 
employ this as a basis to interpret the observed fractionation 
and rotational excitation in a self-consistent fashion (interpret- 
ing brightness, etc. within the context of a model which also 
accounts for the abundance). The following Section sets out in 
detail the physical processes which determine the abundances 
of I^CO and "CO. 



3.2. CO, CH and H2 

By contrast with CO, the variation of N(CH) with N(H2) is very 
nearly linear and with considerably smaller scatter and over- 
all range for the actually-measu red N(CH). Fig . 1 at b ottom is 



an updated version of Fig. 1 of iLiszt & Lucas (l2002h and al 



though it shows a somewhat less-perfect correlation in a larger 
dataset (46 vs 32 sightlines), the derived means < N(CH) > / < 
N(H2) >= 4.1 X lO"** and < N(CH)/N(H2) >= 4.5 + 1.6 x lO"** 



(for lines of sight with CH detections atN(H2) > VxlO'^cm"-^) 
are nearly unchanged, N(CH) and N(CO) are well-correlated 
(coefficient 0.86; this is somewhat contrary to a remark by 
Crenny & Federmaa (2004)) though the functional relationship 
is quite steep, N(CO) cx N(CH)2 <>. 

In general, for N(CH) > 2x 10'^ cm"^, CH is a quite reUable 
indicator of the H2 column density, with N(H2) ~ N(CH)/4 x 
10"^. However, and somewhat remarkably, for N(H2) < 7 x 
10'^ cm"^ it is somewhat more likely to find CO than CH, 
and several lines of sight at lowe r N(H2) have very high valu es 
of N(CO)/N(CH) (see Fig. 18 in ISonnentrucker et al.l (l2007h '). 
Lines of sight with low molecular abundances, including a few 
of those w ith smaller or undete c ted N (CH) in Fig. 1 were dis- 
cussed by IZsargo & FedermanI (l2003h . who argued that non- 
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N(H2) [cm-2] 

Fig. 2. Variatioti of X(i2cO)=N(i2cO)/N(H2) with N(H2). The 
regression line shown has slope 1.014±0.13 and passes through 
XC^CO) = 4.05 X 10-^ forN(H2)= lO^" cm-^ . Lines of sight 
having uv absorption measurements of both '^CO and '^CO 
are noted (red in appropriate media); they have relatively high 
X(i2cO) at a given N(H2). 



thermal processes must preferentially dominate the chemistry 
along some very transparent lines of sight. Alternatively we 
noted that very high values for X(OH) - N(OH)/N(H2) might 
be expected in mo re diff i ise Ho -bearing gas just from quiescent 
thermal processes (ILisztl 120071) : if CO is formed from OH, high 
N(CO)/N(CH) ratios might also be explained in this way. 



4. Mechanisms of CO formation and fractionation 

In this Section we set out the various physical processes con- 
tributing to the formation, destruction and fractionation of car- 
bon monoxide, albeit in a very reductive fashion. The ambient 
gas is taken to be diffuse but purely molecular and the molec- 
ular abundances are taken to small enough that conservation of 
nuclei need not be explictly observed; that is, the fraction of C 
in CO is not large enough to affect the relative abundance of 
HCO^ and the fraction of '^^C in '^CO is not large enough to 
alter the relative rates at which '^CO and '^CO form from the 
isotopic variants of HCO^. 

4.1. Source function for carbon monoxide formation 

As noted above, a source function for carbon monoxide may 
be approximated as the quiescent thermal recombination rate 
of a fixed relative abundance X(HCO^) = 2 x 10"^ such that 
the observations are explained at densities typical of diffuse 
clouds, even if the exact formation route for HCO^ in diffuse 
gas is problematical. Ion-molecule reactions in such a quies- 
cent gas form too little HCO^ and/or CO by a factor of about 
30 and to account for the discrepancy it has been suggested 
that chemical reactions with C^ are driven at higher rates for 



various reasons. The carbon isotope exchange reactions also 
involve C^ and it is of interest to ask whether they also must be 
driven at non-thermal rates in order to account for the observed 
N(i2cO)/N(i^^CO) ratios. 

The rate constant for thermal recombination of HCO^ with 
electrons is 3.3 x 10' ^ cm^ s'VT k according to the UMIST 
reaction rate database (IWoodall et al.. 2007); therefore the vol- 
ume formation rate of '^CO due to recombination of H'^CO^ 
in a gas having an electron fraction X(e)= n(e)/n(H2) is 



dnC^CO)/dt = 6.6 x 10""^ X(e)n(H2)^ cm"^ s^'/Tk 

and the formation rate of '^CO due to recombination of 
H'^CO^ is here taken to be 60 times smalle r following our 



isotope measurements in mm-wave absorption (ILucas & Liszt . 



19981) 



4.2. Carbon isotope excliange 

'^CO molecules are interconverted to '^CO with a rate constant 
kf by the reaction '^C^-H '^CO ^ ''C^-H '^CO + 34.8 K, and in 

kfexpr-34.8/T K) dWatson et al 



the other direction at rate ki 

I — ^— I - I 1 

1976). The rate constants were measured by ISmith & Adamsl 



(1980 ) and shown to be strongly temperature-dependent below 
500 K with a measured value kf = 7 x 10"'° cm^ s"' at Tk= 80 
K, (and an implied value 70% larger at 10 K), as compared with 
kf = 2 X 10"'" cm^ s~' nieasure d at Tk= 300 K bv AVatson etal] 
(ll976l) . ISmith & Adamsl (Il980t) discuss the means by which to 
convert their measurements to values of kf below 80 K and their 
results have been employed b y most subsequent auth ors, for in- 



Chu&Watsoiil(ll983h or lLanger et akl (Il984 (although 



stance . , 

not by Sheffer et al. ~('2007') who used the much smaller value 
of Watson et al. (.1976.) . claiming that it agreed better with ob- 
servation). 

A suita ble expression for kf(TK) based on the results of 



Smith & Adams (1980.) does not exist in the literature. For fu- 



ture reference, we provide the following: 



kf = 7.64 X lO"** Tk"°-^^ cm^ s"' 



(Tk = 80 - 500K) (4a) 



1.39x10-'* Tk"°°^ cm^ s"' 



(Tk = 10 - 80K) (4b) 



1 -Hexp(-34.8/TK) 
In terms of volume formation rates, the forward reaction 
forming '^CO proceeds at a volume rate 

dn(i^CO)/dt = kf XC^CO) X('^C+) n(H2)^ cm"^ s"' 

It is straightforward to show that the backward reaction 
between '^C^ and '^CO is never an important "source" of 
'^CO (given the rate at which '^CO putatively forms from 
H'^CO"^) but conversion of '^CO to '^CO is the dominant 
route to '■'CO formation when the relative abundance X('^CO) 
substantially exceeds 10"^. The '^C-insertion reaction would 
dominate at much smaller X('^CO) if the direct formation 
rates were slower, leading to problems m aking mode l s with 
X('2C0)/X('3C0) » 60, as in Fig. 16 of iWarin et all (Il996l) 
(at far left). 
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Fig. 3. Left: Variation of NC^COVNC^CO) with N (H2) (left) and N('^ CO) . Where possible poi nts are labe led with estimates of 
Tk fr om H2 (at left) and C2 (right) as tabulated bv^ Burgh et a l.' (2007) and Sheff er et al. (2007) for H2 and Sonnentrucker et al 



( l2007i) for H2 and C2. Radiofrequency data are those of iLiszt & Lucas (1998). There is no tendency for the column density ratio 
to vary monotonically with the H2 temperature and only a very loose trend with that of C2 



4.3. Destruction of carbon monoxide by He* 



4.5. Relative abundances of ^^CO and "CO 



The reaction He^ + CO — > + O + He proceeds with a 
rate constant kne^ 1.6 x 10"^ cm-* s"' and is the dominant 
chemical destruction mechanism for '^CO whenever n(H2) 
< 100 cm same is not true for '^CO, whose chemical 

destruction is generally dominated by '^C-insertion whenever 
n(H2) > 2 cm"^, see Eq. 4c and 4d. Calculation of the ion- 
ization equilibrium is complicated, but diffuse cloud models 
tend to produce a nearly constant number density of He^ be- 
cause such a large fraction of the free electrons arises from 
the near-complete photoionization of carbon. From our mod- 
els we take n(He^) - 3.4 x lO"'* cm""* which is appropriate 
either when neutralization by small grains is considered and 
= 2 X 10 s ' per H-nucleon or wh en neutraliza tion by small 



10 '^ s-i ( Lisztll2003h . 



grains is ignored and 

The volume destruction rate of ^^CO by 



He+ 



dn('2C0)/dt = -kHe n(He)+n('2cO) cm^^ s^^ or, numerically 



dn(i2C0)/dt = -5.4 x 10"^^ XC^CO) n(H2) cm^^ s" 



4.4. Photodestruction and selective photodissociation 

The nominal free-space photodissociation rate in the mean in- 
terstellar radiation field is usually taken to be the same for 
either version of CO and we parameterize the photodestruc- 
tion rate internal to a cloud as 2.1 x 10"'" s"' gy I where I 
<> 1 represents the possibility of a variable external radiation 
field, 2.1 X 10"'" s"' is the photodissociation rate in free space 



dLe Teuff et al. ■ I2OOOI) , and gi2 or gn represent the diminution 



of the photodestruction rate due to shielding by dust, H2 and 
other carbon monoxide molecules for either isotope. The free- 
space rates for '^CO and '^CO are take n equal but gi2 declines 
more rapidly into a c loud than does gn dVan Dishoeck & Black 



Given the preceding considerations in this Section we may 
write the following approximate, implicit and local expressions 
for the relative abundances 



X('-^CO) = 



1.2x 10-^X(e)'n(H2)/TK 
gi2l + 0.0026 + 2.3 X 10-%f'n(H2) 



(4c) 



X('3CO) 



(2.0 X 10-''X(e)7TK + 2.3 x 10-5X('2cO))n(H2) 



Where kf' = kf/lO-** — ^ 



;nl + 0.0026 + 0.0014kf'n(H2)exp(-34.8/TK) 

i4d) 

cm" s ', I measures the strength 
of the incident interstellar radiation field normalized to its 
mean strength (1=1) and X(e)' = X(e)/4 x lO""* ; the elec- 
tron fraction is comprised of a contribution from fully once- 
ioniz ed atomic c arbon at the level of 3.2 x 10"^ for an H2 
gas dSofia et al.l I2OO4I). and the re mainder from due to 
cosmic-ray ionization dLiszti l2003l) . The first term in the nu- 
merator of Eq. 4d coiTesponds to recombination of H'^CO^ 
with X(H'2C0+)/X(Hi3C0+) = 60 and the numerical constant 
in both denominators represents destruction by interaction with 
He^. In considering these expressions note that the hydrogen in 
diff'use gas is not wholly molecular, probably even when CO is 
detected, and the total density n(H) - n(H I) H- 2 n(H2) - 2 
n(H2)/fjj2 where fjj^ is the fraction of H-nuclei in H2. 

Note the following with regard to these expressions. For 
'^CO, photodestruction is dominant until gi2 I ~ 0.003, at 
which point the shielding would be strong enough to permit 
near-complete conversion of carbon to CO. Thus this simple 
chemical scheme can caiTy the gas from the diffuse to the dark 
regime (where the free carbon abundance is somewhat lower 
and X('2C0) ~ 10 '*). Also for '^CO, destruction by conver- 



(fl988l) : IWarin etai] ([199©). In Sect. 5 we show an example ^° ^O dominates over destruction by He^ only for n(H2) 



where this behaviour is parametrized as go = g^^f- 



> 100 cm"^; this stands in opposition to the situation for '-^CO, 
where the interaction with He+ is negligible at almost all den- 
sities. This imbalance contributes to the lack of equilibration of 



6 



H. S. Liszt: Formation, fractionation and excitation of carbon monoxide in diffuse clouds 



150 

O 
O 
CO 

^100 
O 
O 

CM 



50 




1=1/3 



^radio 



10 



-6 



X(12C0) 



Fig. 4. Variation of N(i2cO)/N(i3CO) with X(CO) = 
N(^^CO)/N(H2) compared with the model chemistry (curves) 
described in Section 4.5 of the text. Models are for 1=3, 
Tk=45 K, n(H2)= 50 cm^^ 1=1, Tk=45 K, n(H2)=125 
cm-^ 1=1/3, Tk=22.5 K, n(H2)=125 cm^^ (upper of two 
curves) and 250 cm"-*. As in Fig. 3 optical data are labelled 
with the kinetic temperatii re derived from C2 as tabulated by 
Sonnentrucker et alJ OOOtI) . 



the carbon isotope exchange, complicating the interpretation of 
the observed '^CO/'^^CO ratios. 

Finally, note that conversion from '^CO becomes the dom- 
inant source of '^CO only when XC^CO) > 10 '♦/Tr or 
X('^CO) > 2 X 10"^. If the direct formation rate of carbon 
monoxide were taken to be much smaller than that given here, 
as is the case in models which substantially fail to reproduce 
the overall carbon monoxide abundance, conversion from '^CO 
would dominate the formation of '^CO at far smaller X('^CO). 
This would make it very make it very difficult to reproduce 
ratios N('^CO)/N('^CO) » 60, which require strong self- 
shielding of the '^CO and therefore, high X('^CO) (see Sect. 
5.2). 

5. Observed fiC^CO) and fiC^CO) 

5.1. X(CO) at the extremes of abundance 

Figure 2 shows the observed variation of X('^CO) with N(H2); 
with much scatter, t he regression line ha s unit s lope (see the 
caption); Fig. 20 of ISonnentrucker et al.l (l2007h appeared to 
show that X('^CO) » 5 x 10"^, independent of N(H2), for most 
lines of sight at N(H2) < 3 x 10^" cm"^. 



5.1 .1 . The low-abundance limit 

Considering the smallest observed CO column densities and 
abundances, at N(H2) < 3 x 10''* cm^^, <X(i2cO)> * 8 x 
10"** (Fig. 2). In light of Eq. 4c, we note that n(H)/TK ~ 1 in 



diffuse gas, for instance, n(H) = 50 cm"-*, Tk= 50 K yields a 
typical thermal pressure in the ISM, but n(H2)/n(H) < 1 in such 
a gas. Thus even our high ad hoc formation rate for CO predicts 
X('^CO) somewhat below 10"^ in unshielded regions gi2 = 1 
and the observed CO is likely to be largely unshielded at the 
lowest observed N(H2). Gas observed at N(H2) > 10^" cm~^, 
X('^CO) > 10"^ must generally be substantially self-shielded 
in '^CO, implying that the observed N('2cO)/N('^CO do not 
directly reflect the interstellar isotope ratio. 



5.1.2. HighX(CO) 

Complete conversion of carbon to CO in a fully molecular 
diffuse gas would yield XC^CO) = 3.2 x 10"* given the 



free ga s phase abundance of carbon determined bv lSofia et al 



(120041) . implying that the highest observed fractions of carbon 
in carbon monoxide in the molecular portion of the gas are 
['^CO]/['^C] X 0.08. These fractions are small enough to en- 
sure that the observed gas may truly be considered diffuse, with 
nearly all carbon in the form of C^, but they are not necessarily 
small enough that their effects on other species are completely 
ignorable. 

For '^C, the fraction in '^CO may be substantially higher 
owing to fractionation, reaching 0.20 - 0.25 (see just below). 
In this case, the ambient molecular gas is appreciably deprived 
of free '^C-nuclei, and account should be taken of such effects 
on the isotope ratios derived from other species like CH which 
are expected to share the same volume as CO. Such sharing is 
somewhat less obvious for CH^ but in general, it seems clear 
that carbon isotope measurements in the diffuse ISM should 
avoid sightlines having substantial CO, unless the bias caused 
by carbon monoxide fractionation can be corrected somehow. It 
seems generally to be the case that the more pronounced effect 
of CO fractionation is to deprive the ambient gas of '^C but, if 
it were observed that N('^CO)/N('^CO) » 60 at sufficiently 
high X('^CO), the opposite would be the case. 

Again considering Eq. 4c, but in the limit of large X(CO), 
X(CO) ^ 4.7 X 10 VTk in the Hmit of high n(H2) and strong 
shielding ( I gi2 ^ 0) or XC^CO) ^ lO"* at Tk = 40 K. Noting 
the behaviour in Fig. 1 and 2 we infer that this high-abundance 
Hmit would occur only for N(H2) » 3 x 10^' cm"^, N('^CO) =s 
3 X 10'^ cm"^, wefl above the observed range. X('^CO) = lO""* 
is still a factor of three below the free-carbon fraction relative to 
H2 in diffuse molecular gas. In the observed diffuse sightlines, 
the destruction of '^CO is dominated by photodissociation so 
that Igi2 ^ 0.01 even at the largest N(H2) and N('2cO). 



5.2. Fractionation 

Figure 3 shows a summary of available measurements of the 
two most abundant forms of carbon monoxide observed in ab- 



sorption in diffuse clouds. ISheffer et al.l (120071) noted that lines 
of sight studied in '■'CO have about 20 K lower H2 J=l-0 ro- 
tational temperatures than the mean for all H2 surveyed. As 
indicated in Fig. 2, the sightlines examined in '^CO in uv ab- 
sorption are nearly always those with the very highest X('^CO) 
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Fig. 5. Variation of '-CO J=l-0 rotational excitation temperature plotted against N(H2) (left) and N('^ CO). Radio data are 
sho wn as lower limits because they involve the assumption of a beam efficiency. Data from the work of iBurgh et alJ (l2007h 
and lSonnentrucker et alj (l2007h are shown separately; sightlines with measurements in both references are show n chained, note 
that th e former typically derives slightly stronger excitation in these cases. Data labelled "radio" are those of iLiszt & Lucas 
(119981) . Superposed on the data are calculated results for uniform-density spherical models having total densites n(H) = 64 and 



128 cm"^ hke those used to determine N(CO) in Fig. 1, where only a handful of the sightlines represented here require n(H) 
< 64 cm"^. The lower curves at each density employ an older but likely more nearly correct set of rotational excitation rates for 
H-atom + CO collisions (see Sect. 6 of the text). A minimum excitation temperature of 4.0 K is required to produce a /I2.6mm 
emission line having a brightness temperature Tb= 1 K. 



at a given N(H2), and to that extent are not entirely representa- 
tive even though they vary widely in column density. 

The data show N('2cO)/N('^CO) ratios in the range 15 < 
N('2cO)/N('3CO) < 170 and a tendency for the ratio to de- 
cline with increasing column density and/or relative abundance 
X(CO), but with much scatter There is no tendency for the ra- 
tio to decline with the rotational temperature of H2 as noted 
in the Figure but the two lines of sight with the smallest ratios 
studied optically are those with very much the lowest tempera- 
ture indicated by C2. The radio data consistently find smaller 
ratios at a given column density but the effect is more pro- 
nounced in the panel at right. There has been something of a 
recent convergence between the radio and uv absorption stud- 
ies in that the latter have only now found lines of sight where 
N('^CO)/N('^CO) « 60, as was often found to the case in 
mm-wavelength absorption dLiszt & Lucaslll998h . Previously, 
optical absorption studies occured only along lines of sig ht hav- 
ingjnuch larger ratios (Lambert et al., 1994; Federman et al. . 
I2OO3I) and this disparity was the source of some concern. 

The expected functional behaviour of the column density 
ratio is as follows; at very low N(H2) and/or N(CO) the gas 
is unshielded, the photodissociation rates are assumed equal 
and the observed N('^CO)/N('-'CO) ratios should be very near 
the ratio of formatio n rates, presumably ['^C]/['^C] = 60 



dLucas & Liszll Il998l) . However, as indicated in Fig. 2, the 



lines of sight chosen for study of the isotope ratio have such 
high CO abundances that they do not sample the unshielded 
regime. At intermediate X(CO) and N(CO), self- shielding in- 
creases the contribution of isotope exchange to '^CO formation 
at X('^CO) > 2 X 10"^, but also increases the disparity between 
gi2 and gi3 . Therefore, at intermediate X(CO) the column den- 



sity ratio N('^CO)/N('-^CO) can be expected to be both above 
and below the intrinsic isotope ratio in the gas, with the very 
lowest values at larger N(CO). Finally, at very high N(CO), the 
column density ratio must tend toward the intrinsic isotope ra- 
tio when all '^C resides in '^CO but this regime again is well 
beyond the scope of the present dataset. 

In order that the forward and backward isotope insertion 
reactions equilibrate, the rightmost terms in the numerator and 
denominator of Eq. 4d must dominate. In the numerator of Eq. 
4d this implies X('^CO) well above 2 x lO"*" or that gi2 < 0.1 
in Eq. 4c; photodissociation remains the dominant mechanism 
of CO destruction well into the regime where fractionation is 
important. For gnl = 1/33, n(H2) = 40 cm-^TK = 40 K, 
X('^CO) = 4x 10 If the disparity between the isotopic shield- 
ing terms is a factor of a few, Igi3 « 0.07 - 0.1, and dominance 
of the righthand term in the denominator would require n(H2) 
>> 150 cm"^. Thus, although there is substantial creation of 
'^CO via isotope exchange , it does not occur in a portion of pa- 
rameter space where the isotope insertion actually equilibrates, 
creating a reliable gas thermometer 

5.3. What do we learn from ^^CO/"CO ratios? 

Figure 4 shows the observed CO abundance ratios plotted 
against X('^CO) along with the results of some toy models fol- 
lowing the expressions in Eq. 4c and 4d, along with the addi- 
tional parametrization go = g^^^ . The parameter having the 
greatest influence on the observed isotope ratio is the external 
photodissociation rate, rather than the density or temperature. 
In detail, these observations are very hard to interpret in terms 
of either the detailed physical properties of the gas or the in- 
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trinsic carbon isotope ratio. The old hope that CO fractionation 
would serve as a reliable thermometer is not fulfilled. 

However, in very general terms it does appear that the iso- 
topic abundance ratios seen in carbon monoxide can be un- 
derstood in terms of relatively mundane chemical and photo- 
processes in quiescent gas of moderate density, even if it is 
difficult to derive the underlying physical conditions from the 
ratios themselves. This is analogous to the recognition that CO 
can form in the observed quantities from the quiescent thermal 
gas phase electron recombination of the observed amounts of 
HCO^, and in very great contrast to considerations of the abun- 
dances of HCO^ and most other molecular species observed in 
diffuse gas, certainly all the polyatomics and many diatomics 
like CS, whereby quiescent models fail by very large factors. 



6. J=1-0 rotational excitation and line brightness 

6.1. Observed rotational excitation 

In optical absorption, carbon monoxide rotational excitation 
temperatures are derived directly from measurer nents of col 



Sonnentrucker et al 



umn densities in individu a l rot ational levels dBurgh et al. 
20071 



2007i) and N(H2) is generally 
known. At mm-wavelengths, the measurements are of op- 
tical depths and line brightnesses in rotational transitions 
dLiszt & Lucasl 1 1998b and the H2 column density is inferred 
indirectly, by assuming N(HCO+)/N(H2) = 2 x 10"^. 

Shown in Fig. 5 are J=l-0 rotational excitation temper- 
atures Tex(l-O) from these references. They are very nearly 
unanimous in showing (at left) that lines of sight with N(H2) 
^ 5 X 10^° cm-2 have Tex(l-O) < 4 K; 4.0 K is the smallest 
excitation temperature capable of producing a 1 .0 K brightness 
temperature above the cosmic microwave background (CMB) 
in the /I2.6mm J=l-0 line. The same data are plotted against 
N('^CO) at right in Fig. 5, where some disparity occurs be- 
tween the radio and optical absorption data; a few lines of sight 



having N( CO) ~ 10 cm exhibit somewhat brighter J= 1-0 
lines than would be allowed by a 4 K excitation temperature. 

Figure 6 shows the observations recast as integrated 
/I2.6mm J=l-0 rotational brightness temperatures Wco = 
JTb(1 - 0)dv. At radio wavelengths, these are directly ob- 
served to within a scale factor, the beam efficiency. For optical 
absorption data, the line brightness can be crafted from the rota- 
tional level populations, which fix the integrated optical depth; 
the central optical depth is then specified by the b- parameter, 
so that the line brightness can be integrated over the profile. 
For optically thin lines the integrated brightness is independent 
of the b-value, but many of the lines of sight with N('^CO) 
> 10'^ cm"^ observed in uv absorption are predicted to be 
somewhat opaq ue at /I2.6mm, in keep ing with direct measure- 
ments shown bv lLiszt & Lucasl (Il998l) . For Tex(l-O) = 2.73 K 
or 4.0 K, the integrated optical depth of the J=l-0 transition is 
1 ki ns ' whenN('^CO ) = 1.0 or 1.6 x 10'^ cm 2. Finally, note 
that lBurghetal] (l2007h pubhshed only N('^CO) and Tex(l-O); 
to specify the rotational level populations fully we assumed 
that the excitation temperatures of higher-lying levels were the 
same as that of the J=l-0 line. The excitation is weak enough 



that this is not a major problem, little population exists in levels 
above J=2. 

At left. Fig. 6 shows the integrated brightness tempera- 
ture and H2 column density for the optical absorption data. 
Also show n are t he pre dicted d Sonnentrucker et aL , 2007 ) and 
observed ( iLiszt , Il997h brightnes ses toward ( Oph; t he ex 



citation temperature derived by [Sonnentrucker et al.l ( 2007 ) 



is noticeably s maller than earlier values ( Smit h et al 



1978 



Wannier et al.L Il982h leading to a low predicated Wco and 
larger implied optical depth in the J=l-0 line. Corresponding 
data comparing WC O and N(HCO^) were not provided by 
Lucas & Lisztl (119981) . so, to test the consistency of the opti- 
cal and mm-wave results, as elsewhere in this work, we used 
a set of unpublished, very sensitive, low galactic latitude CO 
emission and HCO^ absorption profiles and simply compared 
the integrated CO brightness and HCO^ optical depths, scaling 
the latter ap propriately to form N(H CO^) and N(H2) as done by 
Lucas & Liszt ( 1998) and lLucas & Liszt ( 1996). This is the one 
exception noted in Sect. 2 to the use of previously-published 
data. In any case, it can be seen from Fig. 6 at left that the radio 
and optical data show the same behaviour. 

This comparison of Wco and N(H2) is tantamount to speci- 
fying the CO-H2 conversion factor; indeed, this regime is really 
the only one in which the conversion factor may be derived 
from direct measurement of the two constituents. Clearly the 
ratio of N(H2)/Wco is typically very large for weak-lined dif- 
fuse CO-bearing gas, compared with the usually assumed local 
values N(H2)AVco = 2 - 3 x 102"H2/K-kms-i. However, for 
sightlines of sight with Wco > 1 K kms"\ the ratio is nearly 
canonical or perhaps a bit low. This insight, that the CO-H2 
conversion factor attains nearly its canonical values at rather 
small N(CO) was at the heart of the original discussion of lLisztl 
(1982). 

The plot of WCO vs. NC^CO) at right in Fig. 6 shows 
a well-defined propo rtionality, as noted earlier (Fig. 12 of 
Liszt & Lucasl (Il998l) ) and there is good agreement between 
the optical and radio data when CO emission was actually de- 
tected at A 2.6mm, at N('^CO) > 6 x lO''* cm"^. Apparendy, 
Wco is a fairly robust and accurate estimator of the CO col- 
umn density itself and scatter in the CO-H2 conversion factor 
arises from the vagaries of the chemistry. Note that Fig. 6 at 
right reverts to use of the previously published mm-wave data. 

ForN('^CO) * 3 - 5 X 10'^ cm"- the brightness saturates at 
Wco ~ 3 - 6 K km s ' as the lines become very optically thick. 
This highlights the fact that the transition to complete C^CO 
conversion must occur over a relatively small range of N(H2); if 
diffuse cloud lines of sight commonly show J=l-0 lines of sev- 
eral K, not grossly different from typical dark gas even when 
only a few percent of the free carbon is in carbon monoxide, 
and if the CO-H2 conversion factor is about the same for dark 
and brighter-lined diffuse gas, it follows that there must be 
a very rapid increase of X(CO) over a very small interval in 
N(H2) and Wco, over and above that seen in Fig. 1 . 
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Fig. 6. Integrated '^CO J= 1 -0 /I2.6mm brightness temperatures predicted from uv absorption and/or observed at /I2.6mm, plotted 
against N(H2) (left) and N('^CO) (right). The shaded line at left represents a CO-H2 conversion factor N(H2) = 2 x 10^° cm"^/K- 
kms~'. Dire ct radio ob servations toward 4" Oph are also shown chained to the equivalent optical datapoint. Data fr om the work of 



Burgh et alj (12007.) and lSonnentrucker et alj (120071) are shown separately and the radio observations are those of dLiszt & Lucas , 
1998h at right: see Sect. 6.1 for an explanation of the radio data plotted at left. 



6.2. Interpretation of the rotational excitation 

Until recently, it seemed appropriate to neglect the contribution 
of excitation by H-atoms, for wh ich the cross-sections were 
calculated to be relatively small dGreen & T haddeus, 1976), 
to assume a largely molecular host gas, and to use the CO 
J=l-0 rotational excitation temperature as a probe of the am- 
bient thermal pressure; it could be shown from calculations 
that there wa s a nearly linear inc r ease of Tex( 1-0) with n(H2)TK 
dSmith et all 1 19781 iLisztl \l9T^ iLiszt & Lucaslll99 8') as long 
as Tex << Tk. This variation of Tex(l-O) with the ambient 
thermal pressure of H2 presumably occurs as a result of the 
small permanent dipole moment of CO since the excitation of 
other species with higher dipole moments is typically sensi- 
tive mostly to density. However, a recent calculation of the ro- 
tational excitation rates of CO by H-atoms yielded the result 
that the per-particle excitation by hydrog en ato ms was much 
stronger than by H2 dBalakrishnan et all I2OO2I) . In this case 
even a small admixture of H-atoms can have a marked effect 
on the excitation of CO (ILisztl I2OO6I) and it is not allowable to 
derive the ambient thermal partial pressure of H2 directly from 
Tex(l-O) by assuming a purely molecular host gas. 

As this manuscript was bei ng revised, it was le arned (P. 
Stancil, privat e communication, Shepler et alj (120071) ) that the 
calculation of iBalakrishnan et al.l ( 120021) had been reconsid- 
ered and that the older, smaller H-atom excitation rates were 
in fact more nearly appropriate. Results for both sets of rates 
are shown here. The choice of excitation rates has several in- 
teresting consequences for the interpretation, as noted below. 

Model results for the CO excitation are superposed In Fig. 
5 and 6. These were calculated by running models ( like those 
used to calculate N(CO) in Fig. 1) at fixed n(H) (as indicated in 
the Figures) and N(H) toward the center, plotting results from 



sightlines having impact parameters ranging from the outer 
edge to center For each model of given n(H) and N(H), the ro- 
tational level populations were integrated for each sightline and 
the result plotted as a point in the Figure. Results for n(H) be- 
low 64 cm"^ are not shown because only a handful of the sight- 
lines represented in Fig. 5 are compatible with n(H) < 64 cm""* 
in Fig. 1 (where they are shown outlined and in green to make 
just this point). 

Shown in Fig. 5 are results for models having n(H) - 64 
and 128 H-nuclei cm"^ using the H-atom exc itation rates 



oflo reen &Thaddeu sl('l976') (as approximate d by IWarin et af. 
and those of Balakris hnan et al.l(l2002h . The models in- 



clude excitation by atomic H and He, and ortho and para-H2and 
the effects of photon-trapping were calculated in the micro- 
turbulent approximation for a b-parameter of 0.8 kms"'. The 
characteristic shape of the upper three curves in Fig. 5 is deter- 
mined by an increase in the fraction of atomic H toward the left, 
and by resonant photon-trapping at higher optical depth toward 
the right. With strong excitation by atomic hydrogen, calcu- 
lated values of Tex(l-O) increase with n(H) at small N(H2) or 
N(i2cO). Models using the older cross-sections for excitation 
by atomic hydrogen do not exhibit such an increase because the 
overall excitation is weak when the molecular fraction is small. 

Only if the smaller H-atom excitation rates are appropriate 
will the models used to calculate X(CO) also do an acceptable 
job of reproducing the observed excitation at like density. The 
excitation temperatures derived from the uv absorption data are 
a bit lower than before, but mostly there is the problem that 
the fraction of gas which remains in atomic form in the mod- 
els (which self-consistently determine the local H2 density) is 
large enough to have a profound effect on the excitation. Lines 
of sight with the smallest N(H2) and X('^CO) are expected to 
have formed CO in the least purely molecular gas and therefore 
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to show the highest excitat ion temperatures, opposite to what 
is observed , if the rates of IBalakrishnan et alj ( l2002h are em- 
ployed. Except perhaps for the anomalous-seeming radio data 
near N('^CO) = 10'^ cm"^ at the right in Fig. 5, the calculated 
rotational excitation is far too high, mainly due to the admix- 
ture of residual atomic hydrogen in the CO-bearing regions. 

The situation with respect to the integrated brightness 
shown in Fig. 6 is clearer, as the models actually account 
quite well for the variation of Wco with N(H2) at left and 
they only slightly overestimate the run of Wco with N('^CO). 
Furthermore, the calculated values of Wco are not strongly 
dependent on which set of cross-sections for excitation by H 
atoms is used. This implies that it may be possible to calcu- 
late the CO-H2 conversion factor reliably, although it should be 
stressed that the models form CO rather art ificially. Moreover, 



if the larger H-atom excitation rates of IBalakrishnan et al 



(I2OO2I) are employed, the disparity in the ability of the mod- 
els to account for both Tex and Wco implies that the calculated 
optical depths are much lower than those which are actually 
observed at /I2.6mm or infeiTed from the uv absorption data. 

7. Loose ends and unanswered questions 

The comparison of CO with H2 in diffuse clouds can now 
be based on a dataset which far exceeds that for any other 
trace molecule. Over the diffuse regime, the run of N(CO) with 
N(H2) actually seems more consistent and understandable than 
that of CH, which suffers from an unexplained decline in its rel- 
ative abundance at N(H2) < VxlO'^cm"^ (see Fig. 1); at higher 
N(H2), N(CH)/N(H2) ~ 4 X 10"** with relatively small scatter 
Further observations of CH and other species are required in 
the regime of moderate N(H2) and nascent polyatomic chem- 
istry. 

On average, the mean N('^CO)/N('^CO) ratio is near the lo- 
cal interstellar isotope ratio ['^C]/['-'C] = 60, but fractionation 
effects cause the N('2cO)/N('^CO) to vary between 20 and 170 
with lower values at somewhat higher N('^CO). The fraction 
of free gas-phase '^C in '''CO is large enough in some cases 
that other species sharing the same volume must be somewhat 
starved for '^C, artificially biasing the ratio of their '^C- and 
'^C-bearing variants. However, the effects of sharply-varying 
fractionation might perhaps be most pernicious for mm-wave 
emission studies, which typically rely on ratios of brightness 
temperatures in '^CO and '^CO, coupled with the assumption 
of a fixed abundance ratio N('2CO)/N('^CO), to infer '^CO 
optical depths, excitation temperatures and column densities. 
Unrecognizable systematic variations in N('^CO)/N('-'CO) 
could wreak havoc with interpretation of such datasets. 

At the present time it is unclear whether fractionation in 
CO can be regarded as anything more than a nuisance: the ob- 
served N('^CO)/N('-'CO) ratios result from several competing 
influences including direct formation, carbon isotope exchange 
and selective photodissociation, none of which dominates to 
the extent that the local temperature or density can be inferred. 
At least in the uv datasets, '^CO/'-'CO ratios have generally 
been measured along lines of sight having only rather high 
X(CO) and over a relatively modest range of N(H2) (see Fig. 
2). This may have introduced some bias but interpreting the ob- 



servations is sufficiently difficult that further discussion of the 
fractionation may not be rewarding per se. However, study of 
'^CO may be rewarding for other reasons, especially at higher 
N('^CO) and N(H2) as more complete carbon conversion to CO 
occurs. 

The convergence of the CO-H2 conversion factor for dif- 
fuse and dark gas well within the diffuse regime (at N(CO) 
^ 10'^ cm'^ in Fig. 6 at left) suggests that the full conversion 
of carbon to CO must occur over a very narrow range of ex- 
tinction, N(H2), etc. just beyond the diffuse regime. This may 
leave very little room in parameter space for a "translucent" 
regime in which neutral carbon is the dominant carbon-bearing 
species. Given that the CO abundance in dark gas is typically 
found to be X(CO) = 8 - 10 x 10 ^ some 3-4 times smaller 
than 2[C]/[H] ^ 3.2 x lO"'* in diffuse gas, we are left to ask just 
where the transition from diffuse to dark gas actually occurs in 
terms of N(H2) and what is the carbon budget in the transition 
regime. 

In the diffuse regime, the CO-H2 conversion factor is ac- 
tually measured and found to attain values N(H2)/Wco = 2 - 
3 X 10^"H2/K-kms"' along rather thin lines of sight where the 
fraction of free gas-phase carbon in CO is only a few percent. 
These lines of sight also have Wco of a few K-kms"', as is 
typical of darker material. Again we are left to wonder how 
much room in parameter space is actually left for a translucent 
regime. 

The observed CO-H2 conversion factors and CO J=l-0 ro- 
tational line brightnesses are actually well-explained by the 
uniform density models whose results are presented in Figs. 
1, 5 and 6. The most notable possible failure in the interpreta- 
tion here is the mismatch between the too-large excitation tem- 
peratures and too-small optical depths predicted by the mod- 
els for the J=l-0 rotational transition (see Fig. 5 at left) when 
the recent calculation of the H-atom -1- CO excitation rates 



by IBalakrishnan et al.l (l2002l) is employed (see Sect. 6.2). The 
models are uniform and rather diffuse and so leave a substantial 
fraction of H-nuclei in atomic form in some cases. Excitation 
by such H atoms was ig norable using ol der excitation rates 
Green & Thaddeus . 1976 ) but the rates of I Balakrishnan et"ar 



2002h are so large that even a slight amount of residual atomic 



hydrogen would have a profound effect. Although the disagree- 
ment between observed and measured excitation temperatures 
seems to point to the need to sequester CO in regions of nearly 
pure H2, the agreement is actually worst for the most diffuse 
gas with the smallest N(CO) and X(CO). Should such gas re- 
ally be expected to be the most purely molecular ? Perhaps the 
neatest way around this problem lies with the apparent recent 
realization that the older, smaller excitation rates by H atoms 
are actually more nearly correct. In this case the models do a 
good job of reproducing the rotational excitation of CO at the 
same densities at which CO forms, even if the ambient hydro- 
gen is not exclusively molecular 

There is little difference in '^CO J=l-0 line brightness be- 
tween dark clouds with XC^CO) = 10-^ N(i2cO) = 3 - 
10 X 10'^ cm"^ and that from diffuse clouds with X('2CO) 
= 2 X N(i2cO) = 3 - 10 X 10'^ cm^^. Although dark and 
diffuse gas may better be distinguished on the basis of observa- 
tions of '^CO (acknowledging possible effects of fractionation) 
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or C'**0, such data are not always available. Given the preva- 
lence of diffuse gas at larger distances from the center of the 
Galaxy and larger distances from the galactic plane, it seems 
worthwhile to ask whether possible confusion between diffuse 
and dark gas has caused misjudgement of the quantity and char- 
acter of molecular gas in our own or other galaxies. 
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